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ON CERTAIN APPROXIMATE FORMULAE FOR 
CALCULATING THE TRAJECTORIES OF 
SHOT. 

TN the postscript to a paper by Mr. W. I>. Niven, “On the 
Calculation of the Trajectories of Shot,” which is published 
in the Proceedings of the Royal Society, vol. xxvi. pp. 268-287, 
I have given, without demonstration, some convenient and not 
inelegant formulae applicable to a limited arc of a trajectory 
when the resistance is supposed to vary as the «th power of the 
velocity. 

In these formulae, the angle between the chord of the arc and 
the tangent at any point is supposed to be always small. The 
index n is not restricted to integral values, but may take any 
value whatever. 

As the proof of these formulae is not altogether obvious, and 
a similar method of treatment may be found useful in other 
problems, I think it may not be unacceptable to your readers if 
I show here how the formulae may be demonstrated. 

Analysis. 

Investigation of formulae applicable to a small arc of a 
trajectory, when the resistance varies as the «th power of the 
velocity. 

Let x and y denote the horizontal and vertical co-ordinates at 
time t, u the horizontal velocity, and p the angle which the 
direction of motion makes with the horizon at the same time. 

Hence the velocity at time t is u sec 1 p, and we may denote 
the resistance by ku ' l (sec <£)“, where k is constant throughout the 
small arc in question. 

Also let / and q denote the values of u at the beginning and 
end of the arc, a and B the corresponding values of p, g the force 
of gravity, T the time taken to describe the arc, X and Y the 
corresponding total horizontal and vertical motion. 


Making p the independent variable, the fundamental formulae 
are— 

/ , du iu , ‘+ 1 / ...... 

(1) — =-(sec (/>)" + 1 ; 

d<p g 


(2) 


dx 

dtp 


(sec pf; 


(3) F. = " ^-( sec t an <t> i 

dcp g 

(4) -(sec <py. 
d<p g 


From the first of these equations— 

1 dn k. 
u n + 1 dip 


(sec p)“ + 1 ; 


and therefore, by integration between the limits <p = a and 
<P = ft 

I I _ kn /"* 

r r ~ IJ ^ 

Also, we have— 


and 


and we wish to compare the two former of these definite integrals 
with the following known one, viz. : — 


X = I, 

f « 2 (sec <p) 2 d<p ; 

■Lj 


Y =li 

f ?^ 2 {sec 0) 2 tan <p dtp; 



T = -j 

f «(sec <pY 2 d<p; 

£ 


= («- 2) 


and the last with- 




dp 


dp — 


:) f a 

- I u 2 (sec <p) n + 1 d(f >; 
J $ 


I t ,\ f* I du T, k{n - 1} 
-- = (« - 1) -- ---dp = — ; 


r 


u (sec p) u 4 1 d<p . 


qn-x pn-L ' ' J pti n dip g 

This may be done by means of the following lemma, which follows immediately from Taylor’s theorem :— 

Lemma. 

If F (<p) be any function either of (p only, or of p and u , where I and if a and 0 be the limiting values of (p in the integral and 
u is a function of <p given by the above differential equation (x), | y = ^(a + 0), then, putting for a moment p = y + a>, 

r F(p)dp = f 1£o “^F (y + m)da, = f ” { F(y) + F'(y> + F"(y)^ + F*'(y)** + F s *( 7 )^ + &c. }(/« 

J -/» 2 0 24 J 

= (a-0)i F(y) + -(o - jS) e F"(-y) + - 1 - (<* - /3) J F""(-y) + &c. f 

l 24 IQ20 J 


where F'(d>) = £ 2p^, F"(<|)) = —3$, &c„ and F(y), F'( 7 ), 
dip d<p“ 

F“(y), &c., are what F(p), F'(<|>), F"($), &c., become when 7 is 
substituted for p, and the corresponding value of a (u 0 suppose) 
is put for u. 

In what follows, the last of the terms above written, which is 
of the 5th order in (o - B), is neglected, together with all terms 
of the same order of small quantities. 

All the definite integrals with which we are here concerned 
are included in the two forms 


In the first place, we will apply the above formula to the 
case in which F(ij>) is a function of i p only, viz. when F(pj = 
(sec p)" + 1 , 

Hence 

F'($>) = (« + i)(sec <p) u + 1 tan <p ; 

F"(p) = (« + 1)[(« + 1 )(sec <J>)“ + 1 (tan <j>) 2 + (sec <j>)* + s ] 

= (« + l)[« + 2(secij>)“+ 3 - n + i(sec <J>) M + 1 ]; 
and therefore, 


I u l ( sec 1 p)"‘dp, and / u l ( sec p) m tan p dp. 

J is J 0 

J“( sec <t>y‘+Wp = (a - /3)(sec 7)" + 1 ) 1 + ” /3) 2 [« +l(sec 7 ) 2 - n + 1] j, to the 4th order inclusive. 

j3)(sec 7)" + ] | i + 11—(e “ ft 2 [» + 2(sec yf - « + i] |, 


a 

Hence 


P 


I kn, 

= —(a 


which gives q when p is known. 

In the next place, let F(^>) = id(sec p) m . 

Hence 

F'(p ) = £Ei. = In 1 - (sec p) m + mu l {sec p) m tan p 
dp dp 


"(p) = F(p)^E-u n (sec p) n + 1 + OTtan<*J ; 
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and «”(sec p}" + 1 + m tan +1 4- F(4>)P u n 1 —(see p) n + 1 + + i)«"(sec <p)” + 1 tan <p + m {sec 

L g J L g d<p g J 

or 

F"(<J>) = F(^))j^—^-« 2w (sec <\>f n + 2 + 2^^22”(sec <f>)* + 1 tan (p + 772 s (sec $) 2 — w 2 J 

[ hV n . bj ~~l 

—— « 2w ( sec <p) 2n +- + -/-(n + i)u’ ! (sec (p) n + 1 tan <£ + w(sec </>) 2 J 


f A 2 / hi \ 

= F(<£)-[ — (/ + «)« 2 "(sec (f>)- u + 2 + -"-(2772 + 7 / + i)«” (sec <t>Y l + 1 tan d> + 772(772 + l)(secd >) 2 - m 2 } . 
L pa p- J 


Since 


2/2* « 4-1 , 

. — - u ^ (sec ^>)« + 1 , 


this last expression may be put under the form— 

F"(<£) = F(<£) /(/ + n) + l(zm +72 + 1) tan $ + m (m + 1) (sec <£) 3 - m 2 j-. 

Hence, by the above lemma, 

f //sec <j>”d<p = (a ~ fi)¥(y) j 1 + 2V ( a " / 3 ) 2 [~ W + n )( ^ + l(zm + 72 + 1) ( ~r') tan 7 + m(m + 1) (secy ) 2 - nr ~\ \ 

J fi t L \ud<pj0 \ud(pJo J } 

— (a - 6 ) 11 Q (sec 7)^ 1 + -jz (a ~ jB) 2 (as before) j- 

where ( *\ denotes what ~ becomes when & = o, or when y is substituted for <p, and Un for 22, that is— 

\ua<p J 0 ud<j> T 

( du \ k n , %« + 1 

— r = - 2/ 0 (sec 7) + 

\ud<p) 0 £ 0 

The factor z/ Q may be eliminated from this expression, and the expression itself simplified, by means of the formula— 


n — t ,n — i 

q f 


= (« - l) 


0 u 


I du j, k(n — l) r a 7 , „ +1 j. 

-/Ti -r/ I( P ~ - -' ( sec 40 

- t+1 ^<p g J g 


for, putting m = <n + 1 in the above expression, we have 

J 1 * 7 t . \U 4- 1 r. / I 


u l (sec >p) n + 1 d<j>~(a - j 8 )u ti 
' 0 

Hence 


(sec 7)” +1 ( I + At“~‘ 8 ) ? ^ ;+ ®<(^) 0 + 3 / ( w + I )(^) 0 tan 7 + »+ 1 « + 2 ( sec t) = -(«+*) 2 J}' 

j'* (sec - fl « (sec «)" ^ 1 d<p, or J* (sec *)" ^ 4 - - ^) 

= (sec 7)" , — 1 1 i + (a - fff - n - l) (-^) tan 7+ /« - n - 1 m + « + 2 (secy) 2 - m - n -!»+*+*] j. 


It will be noticed that the term involving / ‘* u \ has disappeared by this division. 

V ud<p} 0 

Now make r« = 2, and this formula becomes — 

// (^c <pf dp = (sec 7 )» 

put 1=2, then, from 
1 ,_) (cos 7)* - 1 ! I 


-n - l 72 + 3 


7&(72 - /)\y” 1 p n — 1 

Divide throughout by g, and put 1 = 2 , then, from before, 

X =_ 1 —\ ( —- - - ' ) (cos 7)» - 1 \ 1 - - (a - Pf f 4 ( —s'. ) ‘an 7 + (« + 4) (sec 7)'- - n + 3I l. 

- 2) "l 24 L V ut£(p /o JJ 

Similarly, divide throughout by t i, r , and put 1 = 1 , then— 

a - =r t)<“ s r)*-’ j - »>' [ a (t5y,“ r + <” + 41 <■“ - " + 3 ] I 


1 


k(n - 1 ) \<j 


Lastly, let 
so that 
then 
and 
Hence 

also 

and therefore 


yi«-i/ ' " ( 24 

F(</>) = tt l (sec p) m tan <p = f(<p) tan <p suppose, 

/(< p) = Msec <p) m ; 

F'(<p) = f'(<p) tan p + f(<p) (sec <pf, 

F"(cp) = f(<p) tan <p + 2/'(<p) (sec <pf + 2/(<#>) (sec cp)- tan p. 

f°F(p)dp = (o - &)\F{y) + a 1 ! (« - £) s F"(7)} approximately, 

J » , 

= (a - 0) |/(7) tan 7 + *(a - 0) ; [/"(7) ‘an 7 + 2 /'(?) (sec yf + 2/(7) (sec 7 ) 2 tan 7] j ; 

[y(py<t>=(a-I 3 ) \f(.y) + M<*- P)-f”(y)} approximately ; 

J 0 

jj(p)dp 4- J“f(p)dp = tan 7 + p.Ca - • 8 )’’[y^ ( sec ?) 2 + ( sec Y) 2 ‘an 7 J i 


1 which the term involving/"^) has disappeared. 
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Now, since f{<p) = u l (sec<py x , we have, as before 


and therefore- 
Hence— 


■*■*»"*]* 

Ay ) — / m t an 

f(-y) \ad<p) 0 

j ^ F (<p)d<p 4 - J f{<p)d<p = tan 7 + A(® - 8) 2 (sec + m + 1 tan 7^ > 


and in the particular case where /= 2, and m = 2, we have' 

Y 
X 


= tan 7 + *(# - + 3 tan 7 ] 

= tan {7 + A(« - + 3 tan 7 ] j • 


Hence the angle which the chord of the arc makes with the axis of x is— 

du 


y + At® - #) 2 jj 

Multiplying by the value of X found above, tve have— 


+ 3 tan 7 


]- 


y t suppose. 


Y = - J^)(cos 7)» - 1 1 { tan 7 - *(« -- W j 7 ) 2 - 4(sec 7) 3 ] + tan y[n~ I « + 4 (*ec 7 )* - 6(*c 7 ) a 

-nA « + 3j| } } 5 
or 

'k{pA){j-^p^) (cos 7)"- 1 {tan 7- At® - ^ { ( 7)2 " 4 ” _I ] + tan7 [»- 2W + S !S£ ‘ rr 


Ye.- 


Considering 


- u-l n + 3 

— l — l — and a - $ to be small quantities of the first order, the above expressions 

?>i -l pn- 1> H 1 


]}>' 


for 


— — -, X, Y, and T are true to the fourth order. 

q n p n 

du 't, which occurs as a factor in some of the terms of the third order may be put under a very convenient 


The quantity (*), ’ 


\ud<p/ 

form in the following manner. 

We have, by Taylor’s theoremj 


( du\ 


(*} + 


/d 2 u\ 
\d(b 2 ) 0 2 


+ &c. 


“ = <*•> + (kh 

In this make a = §(a — ( 3 ) and - J(« - 0 ) successively ; therefore 

y = U 0 -U«~ «( 5) 0 + «« - " &C ‘ 

Hence we have to the first order of small quantities- 


and 


and 

and therefore 


P - q 

a — fi 


( du\ 

afy/o 

\{p + q) = u 0 ; 

( d IL \ = — 2 (? . T_ fj — to the first order. 

\ud<pf a (p + ?)(« - 0) 


Making this substitution for/" —\ the expressions for X, Y, and T become- 
\ud<p /0 


X = 


k{n - 2) 


( ? Ti-^> cos7 > ,, - 1 { 


n ~ 1 p-q. 


I - At-• ArA - 0 ) tan 7 - !i —- (a - 0 )" [n + 4 (sec 7 ) 2 - « + 3 ] } ; 

j / + '/ 


=^"- 2 ) (7^ " ?M (C ° S ^ 1 { 

= ^Y)(^-^) (C0S7), '" 1 { i - 


tan 7 - l-EA( a - 0)[« - 2 (sec 7 ) 2 - » -It] - A(« - 0 ) 3 tan 7 [» - 2 n + 5 (sec yf -ft- I « + 3 ]j 
/ + ? 


-—! tsAL (a - 0) tan 7 

6 / + y 


24 


(a - 0) 2 [» + 4 (sec 7 ) 2 - « + 


3]} 5 


.and these values are still true to the fourth order, considering ——— and a - 0 to be small quantities 01 the first order as before. 

p + q 

The angle which the chord of the arc makes with the axis of x becomes, in like manner— 


-which is true to the third order. 


7 ~ 7 + 'A V (« - /8) + i(« ~ 8 ) 2 t an 7 > 
‘ ’ p + S 
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The above expressions for X and Y may be transformed by introducing this angle 7 into them instead of 7, thus— 
(cos y ) 71 - 1 = (cos 7)”*" 1 - (n- 1) (cos 7)«- 2 sin 7 &) +K° “ 0) 2 tan 7 ^J 

= (cos <y)« - 1 11 - i (a - 0 ) tan 7 - 2-li ( a - £) 2 (tan 7) 2 j- 

x = k^r-E )“ Jn-^j < cos { 1 ( sec yf-^ }> 

= X tan 7 = (cos 7)”- 2 sin 7"[ 1 - 2-2 („ - /8) a [« - 2 (sec yf - » - 3] } ! 


Hence we find- 

and 


X = 


k(n — 2 


Y = 


k{n 




■“ 2 sin 7 Q ; 


Q being — 1 - —— ( a _ [« - 2 (sec 7) 2 — « - 3 ]. 

24 

Similarly, if 
we have 


7' = 7 + i 


P ~ L 
p + q 


(a - 0 ) + i (a - 18) 2 tan 7, 


and therefore 


(cos 7')”- 1 = (COS7)"- 1 - (« - 1) (cos 7)'» — 2 sin yY it—— 0 . ( a - 0 ) + tfa - j 8) 2 tan 7 1 ; 

L p + q -I 

= (cos y )’ ,x ~ 1 { 1 - ?Lz±EzJ_ ( a - 0 ) tan 7 - —±( a - 0 f (tan 7 ) 2 } ; 

l 6 p + Q 4 J 


T = 


k(n 


(yd-i - ( cos 1 - > ~ L (« - [« - 2 (sec 7 ) 2 - » - 3] | 


=_ 1 ( J: _L At 

k(n - I) \q tl — 3 p*— 1 ) 


cos y')" — 1 Q, 


where Q has the same value as before. 

Hence the values of X, Y, and T are as stated in my postscript to Mr. Niven’s paper. 

Although the method of finding the expressions for X and T given above, is perhaps the plainest and most straightforward that 
can be taken, the following leads to simpler operations. 

Let f[(p) — u l ( sec 

Then j/{cp)dcji = j u l ( sec = jj u l ~ n - lt ^d<p by equation (I) 


+ const. 


Hence 

Now let 

then 

and 


k(l - n) 

F(</>) = y'(r6)(sec <p)” ! = u^sec cp)"‘+"+ 1 , 

F'(< p) = f'(<p)(sec <p) m + mf[<p){sec (p) “tan <p, 

F"($) = sec 4 >)’ K + 2m/'(cp)(sec <p) m tan tp ]+ mf(<p)[m(sec <p)’f(tan <p ) 2 + (sec (J>)®+ 2 ] 


= /"(<(>)(sec <p) m + 2/«/ r '(<£)(sec 0)“tan $ +^/(<J>)[?« + 1 (sec <]!>)“+ 2 - r«(sec ^>)“]. 

Hence, by the lemma, 

f F (<p)d<p = (o-^){F( 7 ) + - ^) 2 F"(7)} 

J is 

= (a - £) (/(7){sec7)“ + Jf(a - ^) 2 (sec7)’"[/''(7) + 2?«/'(7)tan 7 + 7«/(7 )[s>k + l(sec y ) 2 - ?»]»] j- 
= (« - ^)(sec 7)“ j/ty) + 5 \(b - ^) 2 [/"(7) + 2mf (7) tan 7 + mf{y)[m + 1 (sec 7) 2 - m ] ] | 


But from above 


/6(/ - ») 




_ n) = J 


A<P)d<h 


= (a - j 8 ){/( 7 ) + *(« - l3) 2 /"(7)}. 


Hence, by division, 

j F(i p)d<p^ k Q-£-—^p »-/-:») = (sec 7)“ | 1 + 2 \(k - fi) 2 ^ 2f«-(2'tan 7 + m{m + t(sec 7) 2 - |. 
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It will be noticed that in this division the quantity f"(y) has disappeared. 
Now, from above, 

fW) = ® ? {*ec <p) n +\ 


and therefore 

and 

Hence 


_ jdu .-\-{n -f i)tan cp, 
ud<p r 


m = n 


J F (<i>)d(i> ” U ~ ( l l ~ n ) — ( sec y) m { 1 + 2*4(“ - £) 2 [_ 2 ^ m («7^) 0 tan y + 2m n + i(tan 7) 2 + m(in + 1 (sec y) 2 - Jj 

= (secy) ;,z j 1 + .-g\(a - / 3 ) 2 £ 2lm(^^j tan 7 + rn(m + 2 n + 3)(sec y) 2 - m(m +2n + 2)J j- 

Now- make /« + «+!= 2, 

or m — — (n — i), and we have 

see <f>) 5 ~ ~ (cos y)" - 1 j I ~||| : (b~,0) s JT(*- l)(^^) tan y+ (« - i)(« + 4)(sec yf - (n - i){n + 3) J | 


In this make l ~ 2, and l = 1, successively, and we obtain the 
same expressions for X and T as before. 

The case thus treated is not one of mere curiosity, but is 
practically important. From theoretical considerations, Newton 
concluded that the resistance of the air to the motion of pro¬ 
jectiles is proportional to the square of the velocity, and very 
little progress has been made in the theory of the subject since 
his time. Experiments have shown that the relation between 
the velocity of a projectile and the resistance offered by the air 
to its motion is far from being so simple as that given by the 
theory. The most extensive and accurate series of such experi¬ 
ments which we have are those made by Mr. Bashforth by 
means of his chronograph, which measures with the greatest 
precision the times taken by the same projectile in passing 
over several successive arcs in the course of its flight. In 
a summary of his results for ogival-headed shot, struck with a 
radius of ij diameters, given in Nature (vol. xxxiii. pp. 605, 
606), Mr. Bashforth concludes that the resistance may be ap¬ 
proximately represented by supposing it to vary as one power of 
the velocity when that velocity lies between certain limits, as 
another power when the velocity lies between certain other 
limits, and so on. 

Thus, if v denote the velocity expressed in feet per second, 
d the diameter of the shot in inches, 


and w its weight in pounds, 

and if — = c, 

w 

then, when v lies between 430 f.s. and 850 f.s., 

the resistance is nearly = 61*3 c ( : 

\iooo/ 

when v lies between 850 f.s. and 1040 f.s., 

the resistance is nearly = 74*4 c ^ ; 

\ 1000/ 

when v lies between 1040 f.s. and 1100 f.s., 

( v \ ® 

-J : 

1000/ 

when v lies between 1100 f.s. and 1300 f.s., 

the resistance is nearly = 108’8 c ( ; 

V1000/ 

and lastly, when v lies between 1300 f.s. and 2700 f.s., 
the resistance is nearly = 141 *5 c 


Hence the resistance varies nearly as the square of the velocity 
both when the velocity is less than 850 f.s., and when it is 
greater than 1300 f.s., but the coefficient increases from 61*3 
in the former case, to 141*5 in the latter. Also, the re¬ 
sistance varies nearly as the cube of the velocity, both when v 
lies between 850 f.s. and 1040 f.s., and also when it lies between 
1100 f.s. and 1300 f.s., but the coefficient increases from 74*4 in 


the former to 108 *8 in the latter case. Again, for velocities 
which are nearly equal to that of sound in air, the proportionate 
increase of the resistance is much greater than that of the 
velocity. 

Mr. Bashforth remarks that the points of transition from one 
law of resistance to another, as stated above, are somewhat 
arbitrary, but that, if they were changed a little in either direc¬ 
tion, the practical error would not be large. 

Of course, if we had at our disposal much more numerous 
and still more accurate observations, it would be possible to 
represent the experimental results with any degree of exactness 
that might be desired, by subdividing the observations into a 
larger number of groups, so that the limiting velocities in any 
one group should be closer together, and that the change of the 
index of the power of the velocity in passing from one group to 
the next should be less abrupt. J. C. Adams. 


SOCIETIES AND ACADEMIES. 

London. 

Chemical Society, December 19, 1889.—Dr. W. J. Russell, 
F.R. S., in the chair.—The following papers were read :—Fran- 
gulin, by Prof. T. E. Thorpe, F.R.S., and Mr. H. H. Robinson. 
The authors prepared the glucoside frangulin from the bark of 
the alder buckthorn (Rhamnus frangula ), and find its formula to 
be C 22 H 22 0 9 . On hydrolysis it yields a yellow product, 
which agrees in its properties with emodin, and a sugar which 
has the power of reducing Fehling’s solution, and is not identical 
with dextrose.—Arabinon, the saccharon of arabinose, by Mr. 
C. O’Sullivan, F. R. S. The substance having an optical activity 
“well above [a]/= 140,” obtained by the author by the 
hydrolysis of arabic acid, and described under the name of a- 
arabinose (Chem. Soc. Trans., 1884, 55), yields arabinose on 
hydrolysis, and appears to bear to this carbohydrate a relation 
similar to that which saccharon (cane sugar) bears to dextrose : 
the author therefore terms it arabinon. It has the formula 
C 10 H 18 O 9 , and on hydrolysis gives a yield of arabinose agreeing 
very closely with that required by the equation T H 2 0 

= 2C 5 H 10 O 5 . As yet it has not been obtained in a crystalline 
state ; it has a specific rotatory power of [a] D = 198° *8, and 100 
parts have the same cuprie reducing power as 58‘S parts of 
dextrose.—On the identity of cerebrose and galactose, by Mr. 
H. T. Brown, F.R.S., and Dr. G. H. Morris. The authors 
give the results of an examination of a specimen of cerebrose, 
prepared from phrenosin, which was placed in their hands early 
in 1888 by Dr. Thudiehum, who first isolated and crystallized 
this substance. They show that its specific rotatory power, 
cupric reducing power, and molecular weight as determined by 
Raoult’s method, are identical with those of galactose, thus con¬ 
firming the recent work of Thierfelder, Zeit. Physiol. Chem., 14, 
209) who has proved the sugar produced by the action of acid on 
cerebrin to be identical with galactose. In the discussion which 
followed the reading of the paper, Dr. Thudiehum said that 
phrenosin, C^H^NOg, consisted of the sugar now shown to be 
identical with galactose, C s H 12 O g , of neurostearic acid, C 18 H 36 0 2 t 
an isomeride of stearic acid, fusing at 84°, and of sphingosine, an 
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